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ABSTRACT 

The power spectrum estimator based on the discrete wavelet transform 
(DWT) for 3-dimensional samples has been studied. The DWT estimator for 
higher than 1-dimensional samples provides two types of spectra with respect 
to diagonal and off-diagonal modes, respectively. The two types of modes have 
different spatial invariance, and therefore, the diagonal and off-diagonal DWT 
power spectra are very flexible to deal with configuration-related problems in the 
power spectrum detection. With simulation samples and the mock catalogues of 
the Las Campanas redshift survey (LCRS), we show 1. the shce-hke geometry 
of the LCRS doesn't affect the off-diagonal power spectrum with "slice-like" 
mode; 2. the Poisson sampling with the LCRS selection function doesn't 
cause more than 1-a error in the DWT power spectrum; and 3. the powers 
of the mass or galaxy random velocity fluctuations, which cause the redshift 
distortion, are approximately scale- independent. These results insure that the 
uncertainties of the power spectrum measurement are under control. The scatter 
of the DWT power spectra of the six strips of the LCRS survey is found to be 
rather small. It is less than 1-a of the cosmic variance of mock sample in the 
wavenumber range 0.1 < A; < 2 h Mpc~^. To fit the detected LCRS diagonal 
DWT power spectrum with CDM models, we find that the best-fitting redshift 
distortion parameter P is about the same as that obtained from the Fourier 
power spectrum. The velocity dispersions for models SCDM and ACDM are 
also consistent with other detections with the LCRS. A systematic difference 
between the best-fitting parameters of diagonal and off-diagonal power spectra 
have been significantly measured. This indicates that the off-diagonal power 
spectra are capable of providing information about the power spectrum of the 
galaxy velocity field. 
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1. Introduction 

The power spectrum is one of most important statistical measure to quantify the 
clustering features of large scale mass density distribution traced by galaxies. It directly 
reflects the physical scales of the processes that affect structure formation. 

In paper I (Fang & Feng 2000), the method of measuring galaxy power spectrum with 
the multiresolution analysis of discrete wavelet transformation (DWT) is developed. It has 
been shown that the DWT analysis provides a lossless estimation of the power spectrum. 
The DWT power spectrum estimator is optimized in the sense that the spatial resolution 
automatically adapts to the perturbation wavelength under consideration. Since the DWT 
analysis has strong capability of suppressing the off-diagonal components of the covariance 
of density fluctuations, the DWT estimator is convenient to deal with data sets with a large 
volume and complex geometry. Moreover, the DWT estimator can fully avoid the alias 
effect which appears in usual binning schemes. Because Poisson process possesses diagonal 
covariance in the DWT representation, the Poisson sampling and selection effects on the 
power spectrum are suppressed to a minimum. 

Besides these technical advantages of the computation of power spectrum, the DWT 
space-scale decomposition provides a physical insight into the clustering behavior in phase 
space, i.e. scale plus physical position. Even for second order statistics, i.e. the covariance 
of density fluctuations, the space-scale decomposition can reveal some clustering behaviors, 
such as two-point correlation in phase space, which can not be measured by decomposition 
with modes non-localized in phase space. Moreover, the covariance in DWT representation 
can avoid the difficulty caused by the central limit theorem (Fan & Bardeen 1995.) 

The algorithm developed in Paper I is mainly for 1-D sample. Many algorithms for 
1-D sample can directly be generalized into 3-D. A problem especially for 3-D DWT power 
spectrum estimator is on rotational invariance. 3-D DWT modes is given by the direct 
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product of 1-D DWT modes, and therefore, each of the 3-D DWT modes is described 
by three positive integers (ji, j2) js) corresponding to the scales of the mode in xi, x-2 
and directions, respectively. Thus, there are two types of modes: diagonal mode with 
ii = J2 = jz] off-diagonal mode for which the three numbers (ji, ^2,^3) are not the same. 
The DWT estimator can provide two types of power spectra: 1. diagonal power spectrum 
given by the powers on diagonal modes, 2. off-diagonal power spectrum given by the powers 
on off-diagonal modes. 

Because the two types of modes have different spatial invariance, the diagonal and 
off-diagonal DWT power spectra are very flexible to deal with conflguration-related 
problems in the power spectrum detection. In this paper, we will study both diagonal 
and off-diagonal power spectra with the Las Campanas redshift survey galaxies, and 
relevant simulation samples and mock catalogues. The purpose is to demonstrate how the 
two types of power spectra are applied to estimate the uncertainties of power spectrum 
measurement caused by various anisotropics of the sample. It includes, for instance, the 
irregular geometry of the survey with very low filling factor, the scale-dependence of redshift 
distortion, and the inhomogeneity caused by the selection function. These results insure 
that the DWT estimator is capable of providing a stable measurement of the LCRS galaxy 
power spectrum and the best-fitting parameters of bias and velocity dispersion. 

The paper is organized as follows. In §2, the DWT power spectrum estimator is briefiy 
introduced with the emphasis on the property and algorithm for the off-diagonal DWT 
power spectrum. §3 describes the real and mock samples of galaxies. The applications 
of the diagonal and off-diagonal power spectra to the configuration-related problems are 
presented in §4. The results of the DWT power spectra of the LCRS, both diagonal and 
off-diagonal, and their error estimation are given in §5. The best-fitting values of the bias 
parameters and velocity dispersion are given in §5 as well. Finally, conclusions are stated in 
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§6. 



2. DWT power spectrum estimator 

2.1. The DWT decomposition of galaixy distributions 

If the position measurement is perfectly precise, the observed galaxy number density 
distribution can be written as 

%(^) = J2 ^mS^i^ - y^) (1) 

m=l 

where coordinate x = {xi,X2,xs), Ng is the total number of galaxies, is the position 
of the mth galaxy, Wm is its weight, and the 3-D Dirac 5 function. Without loss of 
generality, we consider the sample which is contained in a spatial box with spatial ranges of 
xi — — Li, X2 — — L2 and x^ — — L3. It is not necessary that the box L1L2L3 is fully 
filled by sample, i.e., allows ng{x) — within some areas inside the box. 

The observed galaxy number density distribution is considered to be a realization of 
a Poisson point process with an intensity n(x) = rig(x)[l + 5(x)], where ng{x.) is selection 
function, and 5(x) the density contrast fluctuation on the underlying matter fleld. The 
goal of power spectrum measurement is to estimate the power spectrum of the density 
fluctuations S{x) — [n(x)/ri(x)] — 1 from the observed realization (5g(x) = [ng(x)/fig(x)] — 1. 

Different power spectrum estimators adopt different decomposition of Sg{x). The DWT 
decomposition is based on the orthogonal and complete set of 3-D wavelet basis {■0j,i(x)}, 
which can be constructed by a direct product of 1-D wavelet basis as 

where ji — 0,1,2.. (i = 1,2,3) and k = 0...2^^~^. The wavelet V'j.iW is localized 
in phase (scale and physical position) space, i.e. it is non-zero mainly in a volume 
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Li/2^ixL2/2^'2xL3/2^'3, and around the position (xi = hLi/2^\x2 = l2L2/V\x^ = 
(Fang & Thews 1998.) The bases eq.(2) are orthonormal with respect both j and 1, i.e. 

where is the Kronecker Delta function. 

The DWT decomposition of [ng{x.) / ng{x)] — 1 is then given by 

V^j,i(x)(ix, (3) 



Since wavelets are admissible, i.e. / '4'j^i{x)dx = 0, eq.(3) becomes 



~1 ^ f ng[^) 

J J ng{x) 



^j,i(x)dx. (4) 



Thus the fluctuations in galaxy number density can be described by the wavelet function 
coefficients (WFC) ej, which is the density fluctuation on scale j = {ji, j2, Ja) localized at 
the position 1 = {h,h, h}- 



2.2. DWT power spectrum estimator 

For a given scale j, there are totally 2-'i+^2+j3 "WPC variables. All second order statistics 
on the scale j can be drawn from the covariance of these WFCs 

Cf = Jjj!'. (5) 

It has been shown in Paper I that the power of the fluctuations on the modes with the scale 
index j can be estimated by 

^J = ^-^J- (6) 

where 

1 1 2-'l-l 2J2-1 2J3-1 

I? = trC"' = V V V fe!f f7) 

j Ojl+j2+j3 j Ojl+j2+j3 ^ L jJ ' ^ I 

^ ^ il=0 ^2=0 ^3=0 
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and 

2J1-1 2-^2-1 2J3-1 , 2 t \ 

Ns = ^-^ y. y. y. r^d^. (s) 

J 2^,+,,+,3 ^A. ^A. J I ; 

The physical meaning of eq.(6) is clear. The term I? is the mean power of j modes measured 
from the observed realization Ug (x) , and the term Nj is the power on j modes due to the 
Poisson noise. For a volume-limited survey, the mean galaxy density fig is independent of 
the redshift. The Poisson noise power is thus simply l/fig. Pj is usually referred to as the 
DWT power spectrum. 

Paper I showed that for 1-D sample in space L the DWT power spectrum Pj is a band 
averaged Fourier spectrum given by 

-1 oo 

Pj-oJ E \i^{n/n'P{n), (9) 

n=— oo 

where P{n) is the Fourier power spectrum with wavenumber k — 27rn/L (for simplicity, we 
refer to n as the wavenumber below), and ■0(n) is the Fourier transform of the basic wavelet 

Mn) = [ ijMe-'^^'^'^dr). (10) 
Jo 

which is orthogonal to the monopolc, i.e. V'(O) = and localized in wavenumber space. For 
the Daubechies 4 wavelet (Daubechies 1992), |'^(n)p has symmetrically distributed peaks 
with respect to n = 0. The first highest peaks are non-zero in two narrow ranges centered 
at n = irip with width Aup (see Fig. 1). Besides the first peak, there are "side lobes" in 
|-?/'(n)p (sec also Fig. 1). However, for the Daubechies 4 the area under the "side lobes" is 
not more than 2% of the first peak. Therefore Pj is good estimation of the band-averaged 
Fourier power spectrum centered at wavenumber 

rij = np2^, (11) 

with the band width as 

An = 2^Anp. (12) 
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For 1-D sample of Lya forests, the DWT power spectrum is found to be smoothly related 
to the Fourier power spectrum by eqs.(ll) and (12) (Feng & Fang 2000.) 

2.3. Diagonal and off-diagonal modes 

It is straightforward to generalize eq.(9) to 3-D. We have 

1 oo oo oo 

Pi-^j^TJ^TTs E E E |^(ni/2^0^(^2/2^'^)^(n3/2^'^)rP(ni,n2,n3). (13) 

ni=— oo n2=— oo n3=— oo 

Because the cosmic density field is isotropic, the Fourier power spectrum P (711,712,^3) is 
dependent only on 



n = ^Jnl + nl + nj. (14) 

Obviously, the DWT power spectrum is invariant with respect to the cyclic permutation of 
index as 

Ph,h,33 — Pi3,3l,h ~ ^32,33,31 (15) 

Considering eq.(ll) and (14), we can formally define a band center wavenumber rij 
corresponding to the 3-D mode j as 



72j = np^(2ii)2 + (2^2)2 + (2^3)2. (16) 

For an isotropic random field, the Fourier modes with the same n [eq.(14)] are statistically 
equivalent. However, the DWT modes with the same rij [eq.(16)] are not statistically 
equivalent, because the DWT modes are not rotationally invariant. A Fourier mode 
^-i{2n/L){nixi+n2X2+n3X3) ^^^^ ^e obtained by a rotation of mode e-^(2T/i'){nii;i+n^x2+n^x3) 

long as n'l + 7if + n'^ = n\+ n\ + n\. However, the DWT modes don't have the same 
property. Generally, one cannot transform a mode (ji,j2,j3) to (ji,j2)J3) by a rotation, 
even when rij ~ riji. Because of different configurations [eq.(2)] between them, the condition 
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rij — rtji generally does not imply 

= ^j'- (17) 
This invariance holds only when (ji, ^2,^3) is a cyclic permutation of {j'n 32-1 j'z) ■ 

With this property, one can define two types of the DWT power spectra: 1. the 
diagonal power spectrum given by Pj on diagonal modes ji = j2 = J3 = i, and 2. 
off-diagonal power spectrum given by other modes. 

Prom eq.(13), the diagonal power spectrum Pj = Pjjj is related to the Fourier power 
spectrum by 

00 00 00 

Pj= H 12 12 ^jini,n2,n3)P{ni,n2,n3), (18) 

ni=— 00 n2=— 00 713=— 00 

where the window function Wj is 

Wj{n,,n2,n3) = i-|^(ni/2^)V'(n2/2^)V'(n3/2^)r. (19) 
with the normalization 

/oo 
Wj{ni,n2,ns)dnidn2dn3 — 1 (20) 
-00 

It should be emphasized that the window function Wj does not contain any free parameters, 
and is completely determined by the basis used for the DWT decomposition. Therefore, it 
is different from the ordinary convolution with filter containing adjustable parameters. 

Eq.(ll) shows that in the scale space the window function Wj is locahzed around 
n\ — n2 = = np2^. Therefore, the diagonal power spectrum Pj is a band-average of 
the isotropic Fourier power spectrum P(n) with the central frequency n = \ping2K As 
an example. Fig. 1 demonstrates the hnear diagonal DWT power spectra in three typical 
cosmological models, i.e. the A cold dark matter, (ACDM), rCDM, and the standard 
CDM (SCDM). They are calculated by eq.(18), in which the P[n) arc taken to be the 
correspondingly linear isotropic Fourier power spectra (Bardeen et al. 1986). The wavelet 
filter |'0(n)p [eq.(lO)] is displayed in a small frame box. 



- 11 - 



For off-diagonal modes, one can also calculate the linear non-diagonal DWT power 
spectrum Pjij2,j3 cq.(13.) However, in this case, -Pjuajs cannot simply be identified as 
a band average of the isotropic Fourier power spectrum P{n) centered at n = rij. 

Nevertheless, rij is useful to cahbrate the physical scale of a given j. Since we use the 
DWT decomposition, it is convenient to label the spatial scale in figures by j. However, 
the spatial scale for a given j actually is dependent on the size of the sample L, i.e. 
Therefore, for samples with different L, the same j do not have the same physical scale. 
To avoid the possible confusion caused by this, we label the figures on the top with the 
wavenumber k — 2Trnj/L. 

3. Samples 

3.1. The LCRS galaxies 

We use the Las Campanas redshift survey (Shectman et al. 1996) to demonstrate 
the DWT power spectrum estimator, including both diagonal and off-diagonal modes. 
The LCRS consists of 26418 redshifts of galaxies, covering over 720 square degree in 
six 1.5° X 80° strips. Three each of the strips are in the north and south galactic caps, 
respectively. The galaxies in the LCRS are selected from a CCD catalog obtained in 
R-band. The luminosity function has been shown to be best-fitted by a Schechter function 
with M* = -20.29 ± 0.02 + 51ogh, a = -0.70 ± 0.05 and (p* = 0.019 ± 0.001 h^Mpc'^ for 
the absolute magnitude of -23.0 < M - 51ogh < -17.5(Lin et al. 1996.) 

We first convert the heliocentric redshifts of galaxies to comoving distances r in a 
Einstein-de Sitter universe. As the survey depth is not very large, the differences among 
the distances given different models are less than 2%. This correction has been taken into 
account in our calculation. Analysis is carried out for the sample with upper cut-off of the 
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recession velocity 51000kms~^, which is approximately the survey depth ~ 450 h~^ Mpc. 
To perform DWT power spectrum analysis, we further enclose each strip into a cubic box 
with size of 500 h~^ Mpc. Accordingly, the spatial size resolved by one- dimensional mode 
on a scale j is 500/2-^ h~^ Mpc. 

A slice-like geometry of the strips fills only a very small fraction of the cubic volume. 
We define a fiUing factor by 

fs-^ (21) 

where Vs is the volume occupied by the survey sample, and L^ox is the length of the box. 
For a given scale j = (ji, ^2, is), there are totally A''j = 2-'i+^2+j3 modes. Since the 

wavelet ipj^i is localized in physical space, the WFCs of modes located in the empty regions 
will be zero. Consequently, for the scale j = (ji, j2, Ja), only fgNj modes contribute to the 
DWT power spectrum. In this case, the summations overing li, I2 and Z3 in eqs.(7) and (8) 
run only for the fgN-^ modes, and the normahzation factor i/2^i+-'2+^'3 is replaced by fsNy 

The effect of the boundary can be estimated with the so-called "influence" cone 
(Pando & Fang 1998.) That is, if the base il'j,i{x) is well locahzed in the interval Ax, then 
the WFCs iji corresponding to position xq will only measure fluctuations in the interval 
[{xq — Ax)/2^ + 1,{xq + Ax)/2-'+^]. If xq is at the boundary, the WFCs iji measures the 
difference between the density contrasts inside and outside the boundary. Thus, for a zero 
padding of the empty regions, the power on the boundary mode would be statistically 
overestimated. Hence, the uncertainty (most likely overestimation) of the power spectrum 
caused by the boundary modes would be of the order of N^/ fgNj, where is the number 
of boundary modes on scale j. One can also pad the empty regions with galaxies created 
by a Poisson sampling according to the same radial selection function as the LCRS. In this 
case, the difference between the density contrasts inside and outside the boundary is erased, 
so contribute little to the DWT power spectrum and will cause small underestimation. For 
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small scales <^ fs^j: the boundary effect is negligible for both zero padding and Poisson 
sampling. 

3.2. N-body simulation and mock LCRS catalogues 

To access the accuracy of recovering power spectrum with the DWT, we construct 
mock LCRS samples. First we produce ten realizations of the N-body simulation for each 
model of the standard cold dark matter model (SCDM), low density flat model (ACDM), 
and a variant of SCDM model (rCDM). The parameters (l^Oj A, F, o-g) are taken to be 
(1.0,0.0,0.5,0.55) for SCDM, (0.3,0.7,0.21,0.85) for ACDM, and (1.0,0.0,0.25,0.55) for 
rCDM. The linear power spectrum is taken from the fitting formula given in Bardeen et al. 
(1986). We use modified AP^M code (Couchman, 1991) to evolve 128^ cold dark matter 
particles in a periodic cube of length 256h~^Mpc. For the SCDM model, the starting 
redshift is taken to be Zi — 15, and for LCDM and rCDM, Zi — 25, respectively. The force 
softening parameter rj in the comoving system decreases with time as 77 oc l/a(i). Its initial 
value is taken to be 77 = 384 h~^ kpc, and the minimum value to be r/min = 128 h~^ kpc, 
corresponding to 15% and 5%, respectively, of the grid size. For the single-step integration 
of time evolution, we use the "leap-frog" scheme, the total number of integration steps 
down to 2; = is 600. 

We first assume an unbiased galaxy distribution relative to underlying dark matter. 
The mock catalogues are then constructed proceeding in the following steps. 1. We locate 
the observer randomly in the simulation box. 2. Since the survey depth and width of 
the LCRS are larger than our simulations, we replicate the slice sample periodically be 
extending to redshift z ~ 0.2. 3. From the space coordinate of the galaxy, work out its 
redshift and galactic coordinates, for which the angular mask of the LCRS is further applied. 
4. According to R-band luminosity function and the average sampling rates respectively 
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for 50-fiber and 112-fiber subset of the observation data, we determine whether the galaxy 
is admitted into the mock sample using Poisson sampling. Similar to the estimation of 
Fourier power spectrum in the LCRS catalogue (Lin et.al., 1996), the uneven sampling 
rate has httle effect on measurements of the DWT Power spectrum. Finally in step 5, 
the photometric catalog is then generated by assigning each galaxy a specific absolute 
magnitude. 

4. Estimating uncertainties with the DWT power spectrum 

4.1. Diagonal power spectrum 

In order to demonstrate the reliability of the DWT power spectrum estimator, we first 
calculate the diagonal DWT power spectrum for the simulation ensembles in the SCDM, 
ACDM, and rCDM models. The sample is in a 256^ h~^ Mpc^ box. The results are 
displayed in Fig. 2. For comparison, Fig. 2 also plots the non-linear DWT power spectrum 
calculated by Eq.(18), in which the Fourier power spectrum is taken from the Peacock & 
Dodds's fitting formulae (1996.) Obviously, the DWT power spectrum estimator is excellent 
to recover the analytic models over a wide wavenumber range resolved by our simulations, 
i.e. 0.1 < A; < 3 h Mpc'^ 

For the magnitude-limited mock samples of the LCRS in real space with size 500^ h~^ 
Mpc^, the estimated diagonal DWT power spectrum are shown by solid squares in Fig. 3. 
The error bars are given by standard 1-a deviations taken over 10 realizations. Comparing 
with the power spectrum calculated by Eq.(18), one finds that the DWT power spectrum 
estimator can successfully recover the non-linear spectra of the mass fields on scales j > 3, 
i.e., /c > 0.1 h Mpc~^, while the power is systematically underestimated for A; < 0.1 h Mpc~^ 

Recall that the Fourier power spectrum also shows a suppression of power spectrum on 
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large scales. It is generally believed that the suppression of the Fourier spectrum is due to 
the following two reasons: 1. the uncertainty in determination of mean density (Peacock & 
Nicholson 1991; Tadros & Efstathiou, 1996); 2. the finite volume of the survey. 

The large scale suppression of the DWT power spectrum have different reasons. For the 
DWT power spectrum estimator, the overall mean density is not needed. The uncertainty 
of the mean density doesn't affect the DWT power spectrum. As a numerical test, we 
computed the DWT power spectrum in an arbitrary selected thin slab in the simulation 
box with the same filling factor as that of the LCRS. We found that the scatter of the 
DWT power spectra from different slabs is much less than the suppression shown in Fig. 
2. Accordingly, the uncertainty of the overall mean density is not responsible for the 
suppression of the DWT power spectrum on large scales. 

The suppression of the DWT power spectrum on large scales is because the diagonal 
mode (j, J, j) cannot match with the slice-hke geometry of the survey. For the LCRS, the 
average thickness in declination is estimated to be ~ 15 h~^ Mpc, i.e. only ~ 3% of the 
side lengths in the other two dimensions. Therefore, the perturbations on scales k < 0.2 
h Mpc~^, or j < 3, in the dimension of declination are poorly sampled. This leads to an 
underestimation of power spectrum with j < 3 in the direction of declination (§3.1). Based 
on above discussions, for mock samples with the LCRS survey geometry, the diagonal 
power spectrum ji = j'2 = Js = 3 are capable of recovering the non-linear power spectrum 
on scales j > 3, but suffers an underestimation on scales j < 3. This geometric effect could 
also be fixed by the off-diagonal power spectrum as demonstrated in next subsection. 
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4.2. Slice-like geometry and off-diagonal power spectrum 

Assuming that ji points along the direction of declination, j'2 the R.A. and js the 
radial or redshift direction, the modes (ji,j2,j3) which can fit with the LCRS' slice-hke 
geometry are 2 < j2, J3 < 7 and 4 < ji < 7. We use the off-diagonal power spectrum Pjjjj 
with a fixed ji > 4. 

Figure 4 illustrates the off-diagonal DWT spectrum P^jj measured from subsets of 
the sample in a O.IL x L x L slab for SCDM and ACDM models respectively, in which 
redshift distortion and selection effect have not been taken into account. Clearly, it shows 
that the off-diagonal DWT power spectrum Pqjj estimator is capable of reproducing the 
non-linear power spectrum on large scales. This power spectrum is in good agreement with 
the theoretical curves on large scales as /c ~ 0.5 h Mpc~^. 

Figure 5 plots the off-diagonal DWT spectra Pajj and Pgjj in real space measured 
for the mock LCRS samples in the ACDM model. It again shows that the systematical 
underestimation of the powers on large scales is significantly erased. The off-diagonal DWT 
spectra Pqjj and P^jj can be used for recovering the non-linear power spectrum of the 
mass field on scales as large as A; ~ 0.15 h Mpc~^. 



4.3. Inhomogeneous effect of selection function 

The purpose of introducing selection function ng{r) is to eliminate the inhomogeneity 
and anisotropy of mean number density of galaxies caused by the sampling of the survey. 
That is, the distribution with uniformly mean number density of galaxies can be recovered 
by weighting a galaxy at radial coordinate xs — r by factor l/ng{r). 

However, the effects of selection function actually are two-fold. On the one hand, it is 
to eliminate the inhomogeneity of the mean number density of the galaxy distribution. On 
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the other hand, it causes an inhomogeneity in the correlations of the power spectrum. The 
later can be seen from the band-band correlation in the Paper I, which showed that, if the 
selection function is slowly varying, the observed auto- correlation of power Pj is 



^ 2h-l 2^2-1 2^3-1 

(-^J-^j) ^ {PiPi)i + 92(71+72+73) X] X] 



(22) 



where (PjPj)/ is the power-power correlation before the Poisson sampling, and therefore, 
cyclic invariant with respect to (ji, J2, Js), i-e. 



The number ^(xs) in eq.(22) is the mean of ng{x3) in cell (Zi, ^2, Iz)- Since ng{x2) depends 
only on X3, the last term in the r.h.s. of eq.(22) violates the permutation invariance. That 
is, the correlations {Pjij2,hPh,j2,h) ^'^^ longer cyclic invariant. In other words, the 
probability distribution function (PDF) of Pj is anisotropic. 

To estimate this effect, we calculate the off-diagonal power spectra P^j^ and P^^^j vs. j 
for ACDM simulation sample in real space, and the selection function is taken to be that 
for the LCRS galaxies, but in plane parallel approximation. The result is shown in Fig. 6. 
The two spectra Pij i and P^^j show some difference, which is not significant on all scales 
2 < j < 7 within 1-a. Therefore, with the DWT estimator, the effect of the anisotropy 
given by the LCRS selection function can be ignored at least for the power measurement 
on scales similar to P4j4 and Puj- 



4.4. Redshift distortion and random velocity field of galeixies 



Usually, the redshift distortion is considered to have two effects on the power spectrum 
measurement: 1. the enhancing of power on large scales due to the linear effect of redshift 
distortion; 2. the suppressing of power on small scales due to the random motions of 
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galaxies inside virialized groups and clusters of galaxies. The general theory of the two 
effects on the DWT power spectrum is given in Appendix. An interesting result is that the 
effect of the galaxy random velocity field v(x) upon the power on mode j is mainly given by 
the power of the pairwise velocity dispersion of the same mode. Therefore, the DWT power 
spectrum can be employed for extracting information of the scale dependence of galaxy 
velocity fluctuations. 

Recall that in the recovery of the real space Fourier power spectrum from redshift 
distorted measurement, it is usually assumed that the peculiar velocity dispersion is 
scale- independent. However, theory and measurement of galaxy pairwise velocity dispersion 
do not indicate a scale-independent, but weakly scale-dependent velocity field. Therefore, it 
is necessary to estimate the uncertainty caused by the assumption of the scale-independence 
of the velocity fluctuations. 

Let's use the plane parallel approximation. In this case, the redshift distortion is 
determined by the a^s-component of the galaxy velocity field, i.e. t'3(x). Subjecting V3(x) to 
a DWT decomposition, we have 



separation on scale j at position 1. Therefore, it is the pairwise velocity on scale j at 




(24) 



The WFCs of the velocity field, eji, describes the radial velocity difference with spatial 



position 1. The redshift distorted DWT power spectrum 
Appendix §A.2) 



is approximately given by (see 




(25) 



where 5] is determined by the geometry of the mode j 




(26) 
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and py is the power of the pairwise velocity dispersion on scale j, i.e. 

P^-&)v (27) 
Therefore, comparing PP with Pj, we are able to determine py. 

As an example, we study the assumption of the scale-independence of galaxy velocity 
fluctuations. One can show that, if the power of pairwise velocity is scale-independent, we 
have 

^LhJsPLhJs - const, for a fixed J3 (28) 

Therefore, the scale-dependence of the velocity fluctuations can be detected by off-diagonal 
DWT power spectrum with a fixed js (redshift direction). Figure 7 plots Paj^a vs. j in 
redshift space for the ACDM simulation samples. It can be seen clearly from Fig. 7 that the 
DWT power spectrum -P4j,4 in redshift space is almost parallel to the theoretical non-hnear 
spectrum in real space. The difference between the redshift distortion suppressions of 
mode (4,2,4) and (4,7,4) is no more than a factor of 2. That is, the suppressing factor 
1/(1-1- Sy) ~ const, and therefore eq.(28) approximately hold in the scale range of spectrum 
P4,j,4. Similar results are obtained for models SCDM and rCDM models. The assumption 
of a scale independent power of the pairwise velocity seems to be reasonable at least for the 
DWT power spectrum of the LCRS sample. 

Figure 7 also presents the off-diagonal power spectra -P4,4j. Obviously, -P4,4j 
significantly differs from ^4^,4 in both shape and amplitude of the power. Comparing to 
-P4,j,4, the amplitude of i-*4,4j is enhanced on large scales and suppressed on small scales. 

The difference between the spectra -P4j,4 and -P4,4j in redshift space is worth while to 
look at. For an isotropic velocity field, Pj^j^Js cychcally invariant with respect to the 
index (ji, j2, js), i-e., 

PiAj = -f4j,4- (29) 
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The difference between -P4,4j and -P4j,4 is simply given by tfie geometric factors 5j, which 
are not cychcally invariant 

S',,,j > 5;,- 4 i>4; (30) 

Therefore, using the difference between i-*4,4,j and -P4j-,4 in redshift space we are able to 
directly measure the power spectrum of galaxy velocity field, P^^^j or P^j^^- 



5. The DWT power spectrum of LCRS galELxies 

5.1. The diagonal DWT power spectrum 

With the preparations given in last sections, we can measure the DWT power spectrum 
of the LCRS galaxies. As discussed in §4.1, the diagonal power spectrum could provide a 
robust measurement of the non-linear power spectrum on scales of A; > 0.1 h Mpc~^. In Fig. 
8, we present the results for the LCRS, where the diagonal DWT power spectra measured 
in the six slices are plotted by scatter symbols as indicated in the figure. The average value 
over the six strips, and over the three north as well as three south strips are also displayed 
by the connected lines, respectively. 

Although the geometry of the six strips are actually quite different, and the fiUing 
factor is only about 0.9%, Fig. 8 indicates that the scatter of the power on scales of 
0.1 < A; < 2 h Mpc~^ is rather small. Comparing with Fig. 3, we find that the scatter of 
the power in Fig. 8 is within the range of l-cr deviations given by the mock LCRS samples. 
Therefore, most of the scatters in Fig. 8 are probably from the cosmic variance. This result 
implies that the DWT power spectrum estimator is insensitive to the geometry of survey 
and filling factor. This feature allows us to directly compare the DWT power spectrum 
measurements from different surveys, or different region in a survey. 
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As indicated in Fig. 8, the DWT power spectrum is capable of probing the fluctuation 
power on scales as small as = 2 h Mpc~^ (j = 7). This is partially because the DWT 
decomposition is free from the "aliasing effect" which may lead to underestimation of power 
on the scale of grid size (Jing 1992, Baugh & Efstathiou, 1994). At the scale j — 6, the 
number of effective modes is fgNj ~ 2500 comparable to galaxy number in each slice, and 
thus the error bars for j > 6 are mainly from the Poisson noise. Even on the smallest scale 
j — 7 the Poisson correction is still much smaller than the detected power. Therefore, 
the diagonal DWT power spectrum from A; ~ 0.1 to 2 h Mpc~^ is quahfled for model 
discrimination. 



5.2. Fitting values of the bias parameter and velocity dispersion 

In order to find the best values of the bias parameter and velocity dispersion, we 
compare the detected LCRS diagonal DWT power spectrum (Fig. 8) with the fitting 
spectrum in redshift space, Pf, given by eq.(13), i.e. 

1 oo oo oo 

^f^^J^TJ^Tj; E E E |V;(ni/2^-^)V'(V2^'^)V'K/2^'^)rP^K,n2,n3), (31) 

ni=— oo n2=— oo n3=— oo 

where P'^(ni, ^2, ris) is the Fourier power spectrum in redshift space. For typical CDM 
models, P'^(ni, ^2, ria) can be taken from the fitting formula given by Peacock and Dodds 
(1994). It is 

P^{k)^h''Pm{k)G{(3,ka,) (32) 

where Pm{k) is the non-linear mass power spectrum in real space, j3 = Q^'^/b is redshift 
distortion parameter, b the scale-independent linear biasing parameter, a^, is 1-D peculiar 
velocity dispersion. G{P, ka^) describes the suppression of redshift distortion. 

Exactly speaking, one cannot simply generalize eq.(13) to eq.(31), because the redshift 
distortion of the DWT power spectrum is different from the Fourier power spectrum. The 
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pairwise velocity suppression of the DWT modes should be calculated based on eq.(25), 
rather than via the Fourier modes. Moreover, cq.(32) is given by averaging over the 
directions of Fourier modes, and therefore, it doesn't match the anisotropic mode (ji, j2, J3) 
of eq.(31). Nevertheless, the fitting of eq.(31) will provide worth information for diagonal 
modes, as the power of the galaxy random velocity field is approximately scale-independent 
(§4.4). 

We treat b and in eq.(31) as free parameters, and fit the power spectrum eq.(31) to 
the detected DWT power spectrum of LCRS. The result is shown in Fig. 9, in which the 1-a 
error bars of the observed DWT power spectrum are given by the scatter over the observed 
6 shces themselves. As have been discussed in §4.1, the DWT power spectrum at large 
scales (j < 3) are underestimated due to the slice-hke geometric effect. We approximately 
compensate this effect in the observed DWT power spectrum by using the fractional 
damping of amplitude calculated from the mock samples in redshift space relative to the 
expected values from the non- linear fitting formula (Eq.(32)) in each model respectively. 
The least square best fitting parameters b and ay for each of models are presented in Table 
1. 

The SCDM and ACDM parameter f3 given in the Table 1 is about the same as other 
measurements, such as P — 0.52 ± 0.13 for the 1.2Jy IRAS survey (Cole et al. 1995); 



Table 1 

Model b P{Pjjj) PiPbjj) CFv km s~^ {Pjjj) (^v km %~^{P^j^ 

SCDM 1.78 ±0.30 0.56 ±0.09 0.50 306 ± 106 397 

rCDM 1.50 ±0.18 0.67 ±0.08 0.59 135 ± 42 217 

ACDM 1.06 ±0.13 0.46 ±0.06 0.41 250 ± 72 344 
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P = 0.47 ± 0.16 (Tadros et al. 1999) and /3 = 0.4llg:Ji (Hamilton et al. 2000) for IRAS 
Point Source Catalog Redshift Survey (PSCz). Actually, our estimated P in each CDM 
model is very close to the linear value UgOo^. This result is consistent with the theory 
of redshift distortion in the DWT presentation, which shows that the linear effect on the 
diagonal DWT power spectrum is about the same as that given by eq.(32) (see. Appendix 
§A.l.) 

The 1-D velocity dispersions given in Table 1 look lower than other measurements in 
the LCRS galaxies. They are Vi2 = 570 ± 80kms^^ (Jing et al. 1998), and 363 ± 13 km s^^ 
(Landy, Szaley and Broadhurst 1998). Yet, the first measurement is based on a priori infall 
model of galaxies. In the second measurement Vi2 is given by the width of an exponential 
distribution, or a Lorentzian distribution in k space. On the other hand, cr„ in eq.(31) is 
assumed to be the variance of a Gaussian distribution of pairwise velocity. To fit a given 
power suppression by Lorentzian and Gaussian distributions, Vi2 will always be larger than 
ay. For instance, if the suppression factor is equal to about 3 as shown by the spectra of 
Fig. 7. fi2/o'i. will be ~ 1.4. Therefore, the values of given by the fitting of the SCDM 
and ACDM basically are consistent with the second measurement. 

5.3. Off-diagonal DWT power spectrum 

As an example we detected the off-diagonal spectrum P^jj of the LCRS galaxies. The 
result is presented in Fig. 10. 

As mentioned in last section, eq.(31) cannot be used to fit with off-diagonal DWT 
power spectrum, as P^{k) is isotropic. One can expect that the fitting parameters of the 
spectrum eqs.(31) and (32) with an off-diagonal power spectrum would be systematically 
different from those of the diagonal power spectrum. This systematic difference actually 
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contains valuable information. The best fitted spectra of P^jj for models SCDM, rCDM and 
ACDM are also shown in Fig. 10. The best- fitting parameters are given in Table 1, which 
clearly show the difference between the parameters of diagonal and off-diagonal spectrum. 
Now we analyze the physical meaning of the difference. 

From Appendix A.l, the linear redshift distortion of spectrum Pj is determined by the 
factor where Sj is a geometric factor given by eq.(A4). The hnear redshift distortion is 
significant on large scale. Moreover, the geometric factor satisfies 

S^jj < Sjjj for j < 5. (33) 

Therefore the spectrum P^jj is less affected by the linear redshift distortion than Pjjj. 
Thus, if fitting both Pjjj and P^jj with the same fitting formula eq.(31), the best- fitting /3 
for P^jj will be lower than those of Pjjj. Table 1 shows the systematic lowering of P{Pbjj) 
than l3{Pjjj). 

The effect of galaxy pairwise velocity on spectrum Pj is determined by the power P^ S'^ 
(Appendix A. 2). If the power of pairwise velocity is approximately scale-independent, we 
have 

That is, the spectrum P^jj is strongly affected by the random velocity field than Pjjj. Thus, 
if fitting both Pjjj and P^jj with the same fitting formula eq.(31), the best-fitting a.^ for 
P^jj will be larger than those of Pjjj. This is shown in Table 1. The significant difference 
between the best-fitting results of Pjjj and P^jj shows that the off-diagonal power spectrum 
has the capability of measuring the galaxy velocity field from the observational data. 
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6. Conclusions 

The DWT power spectrum estimator is studied with the galaxy distributions of the 
Las Campanas redshift survey, and the relevant mock and simulation samples in models of 
the SCDM, rCDM and ACDM. 

The DWT estimator for 3-D samples provide two types of spectra with respect to 
diagonal and off-diagonal modes, respectively. The two types of the DWT local modes have 
different configuration and invariance, and therefore, the diagonal and off-diagonal DWT 
power spectra are flexible for studying direction-dependent properties. 

Using mock and simulation samples, we studied the effect of the irregular geometry 
of the survey, the inhomogeneity caused by the selection function, the scale-dependence of 
the redshift distortion due to the random velocity field of galaxies etc. To estimate the 
uncertainty and errors caused by these effects, the localization of the DWT modes in phase 
space are essential. For instance, the off-diagonal power spectra are capable of measuring 
the statistical features of the random velocity field of galaxies on mode j, because the DWT 
power on mode j is mainly affected by the power of the velocity fiuctuations on mode j. 

The DWT power spectrum from the six shces of the LCRS galaxies are very stable 
over the scale range 0.06 < k < 2. Oh Mpc~^, in spite of the fact that the filling factor of the 
LCRS galaxies is only about 0.9% (in comparison, SDSS's fiUing factor is about 17%.) To 
fit the LCRS DWT power spectrum with CDM models, we find that the redshift distortion 
parameter (3 and the 1-D velocity dispersion are consistent with results of the Fourier 
power spectrum. Therefore, the DWT estimator can give a robust measurement of the 
banded Fourier power spectrum. 

More interesting, the differences between the best-fitting parameters of LCRS' diagonal 
and off-diagonal DWT power spectra are found to be significant. This shows that the 
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redshift space behaviors of the diagonal and off-diagonal DWT power spectra are different. 
With the diagonal and off-diagonal DWT power spectra, we would be able to compare the 
pairwise velocity dispersion on different modes. Therefore, the difference between diagonal 
and off-diagonal power spectra would be valuable to constrain models. 

LLF and YQC acknowledges support from the National Science Foundation of China 
(NSFC) and National Key Basic Research Science Foundation. We thank anonymous 
referee for helpful comments. 



A. Redshift distortion in the DWT representation 

A.l. Linecir redshift distortion 

If the density field is viewed in redshift space, the observed radial position is given by 
the radial velocity consisting of the uniform Hubble flow and the peculiar motion v(x). In 
the linear regime of clustering, the effect of redshift distortion can be described by a linear 
mapping from density contrast 6 (x) in real space to redshift space 

S^{x) = S5(x), (Al) 

where S is the linear redshift distortion operator. In plane-parallel approximation, it is 

^^^ + ^5^"'' ^^^^ 
where (3 is the redshift distortion parameter, and is on the direction of redshift. 

The differential operator in eq.(A2) is quasi-diagonal when decomposed into the DWT 
basis (Farge et al. 1996.) Thus, similar to eqs.(2) and (3), 5^{x) can be decomposed into 

e']c^il + (3Si)e\, (A3) 



-27- 



where the coefficient 5j is given by 

This integral is independent of 1. For diagonal modes j = {3,3,3), we have 

Sjn = 1/3. (A5) 

Using eqs.(5) and (6), the relation between the DWT power spectra in redshift space, 
and real space is given by 

P^^{l + (3SsfPy (A6) 
For diagonal DWT power spectrum, we have 

Pfn = (1 + \l^fPnr (A7) 

The redshift distortion of the Poisson correction are not considered in eqs.(6) and (7). 
The presice theory of the redshift distortion of the Poisson correction requires to calculate 
the selection function by the DWT decomposition. We will not give it here. In a rough 
approximation, the effect of the Poisson correction can be estimated by 

P^^{l-r PS^fl? - A^j = (1 + PS^nPi + Aj) - Ny (A8) 

A. 2. Effect of random velocity field 

In non-linear regime of clustering, the redshift distortion generally is estimated by 
assuming a random velocity field v(x), which describe the velocity dispersion of galaxies. 
The average of v(x) is zero, (v(x))t, = 0, where (...)„ is for the average over ensemble of 
velocity fields. The covariance of the density contrast in redshift space is then given by 

{{S{xi, X2,X3 + vs(x)/Ho)S{x[,x'2, 4 + V3{x')/Ho))y), (A9) 



Thus, the DWT power spectrum in redshift space is 



where 



6^' 



'1 = y 5{xi,X2,X3 + V3{x)/Ho)llji^i{x)dx. 



Using an auxihary variable J, eq.(AlO) can be rewritten as 



Thus, we have 



~s. 
e 



] = J dx5 (^i,X2, 



6 

X3-^-|e 



J=0 



I 



dxd-x.'5 (xi,X2, X3 - i-^j S lx[,X2, 4 - i-^ 



iJvsix)/Ho+iJ'v3{x')/Ho'. 



J,J'=0 



In the lowest order of non-zero correction, we have 



l^^iJv3{^)/Ho+J'vi{^')/Ho^^ ^ 
1 



1 - 



In this case, the redshift distorted power spectrum is 



where is the WFCs of velocity field, i.e. 



^1 



The geometry factor S[ is positive, and is given by 
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When the correction term is large, we have an approximation as 



If the velocity field is homogeneous and isotropic, (|eyP)t; is 1-independent. It is the 
power of velocity fluctuations on scale j, i.e. 

Py = {\^/)v (A19) 



Thus, we have 

' j ~ i + s(py/Hf' 

Therefore, the effect of the galaxy random velocity field upon the power Pj on scales j is 



mainly given by the power of the vs fluctuations on the same scale. 

Similar to eq.(A8), the redshift distortion of the Poisson correction can be estimated by 

= (l - -^fM/)^) (^J + ^J) - ^J- (A21) 
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Fig. 1. — Linear diagonal DWT power spectrum in three typical cosmological models: 
ACDM, rCDM and SCDM. The wavelet filter is displayed in a small frame. 
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Fig. 2. — The non-linear diagonal DWT power spectrum of simulation samples (solid square) 
and the corresponding DWT power spectrum (solid line) calculated by eq.(18), in which the 
non-linear Fourier power spectrum P{ni,n2,ns) are taken from Peacock & Dodds's fitting 
formulae (1994). The error bars are given by 1-a deviations obtained from 10 realizations. 
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Fig. 3. — The diagonal DWT power spectrum of the mock flux limited LCRS catalogue 
in real space for models SCDM (upper panel), rCDM (central panel) and ACDM (lower 
panel). The non-linear DWT power spectra (solid line) are calculated by eq.(18), in which 
the Fourier spectrum P{ni,n2-,n^) is taken from the Peacock-Dodds fitting formulae. The 
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Fig. 4. — The off-diagonal DWT spectrum Pqjj for SCDM and ACDM simulation slice 
samples in real space. Selection effect is not applied. The dimension in x-axis is thin 
{x — O.IL). The sohd line is the spectrum calculated by eq.(13), in which the Fourier 

spectrum P{ni,n2,ns) is taken from the Peacock-Dodds fitting formulae. 
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Fig. 5. — The off-diagonal DWT spectra P^jj and Pajj for ACDM mock LCRS samples in 
real space. The sohd lines are the spectrum calculated by eq.(13). 




Fig. 6. — The off-diagonal power spectra P^j^ and Puj vs. j for ACDM simulation sample 
in real space, but the selection function is taken account by plane parallel approximation. 
For clarity, the error bars of P^j^ are shifted left by a factor of 0.2. 
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Fig. 7. — The off-diagonal power spectra P^j^ and Puj vs. j for ACDM simulation sample 
in redshift space. The redshift distortion is taken account by plane parallel approximation. 
Selection function is not applied. 
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Fig. 8. — The DWT power spectrum measured in the flux-hmited LCRS samples. The 
scatter symbols represent the DWT power spectrum measured in the six slices as indicated 
on the legend. The solid line show the mean power averaged over these six slices. The dash 
and dot-dash line are for the DWT power spectra calculated from the three north and south 
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Fig. 9. — The best-fitting (solid line) of the LCRS diagonal DWT power spectrum for models 
SCDM (upper panel), rCDM (central panel) and ACDM (lower panel) with parameters hsted 
in Table. 1. The observed values (sohd square) has been corrected for shce-like geometry 
effect using the mock samples. The error bars are given by 1-a variance obtained from the 
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Fig. 10. — The best-fitting (solid line) of the LCRS off-diagonal power spectrum P^jj 
for models SCDM (upper panel), rCDM (central panel) and ACDM (lower panel) with 
parameters listed in Table. 1. 



